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Abstract. We consider a perfectly homogeneous, isotropic and spatially flat 
universe which undergoes a sudden phase transition producing topological defects. 
We assume that these defects form a coherent network which scales like the 
background density during the radiation and matter dominated eras, and describe 
them in terms of a few free functions. We carefully model the loss of scaling invariancc 
during the transition from radiation to matter dominated era. We choose a simple set 
for the free functions, compute the microwave background temperature anisotropies 
generated by such a network, and compare our results to previous calculations. 



1. Introduction 

The origin of the small inhomogeneities in the cosmic fluids which eventually evolved 
into the large scale structures that we observe today in the universe is still under 
debate. Two classes of models are at present in competition. The inflationary 
models on one hand jj] explain those inhomogeneities by the amplification of quantum 
fluctuations at the end of an accelerated phase of expansion of the universe || . On 
the other hand topological defects appearing during a phase transition in the early 
universe H can also seed inhomogeneities 0. 

A way to discriminate between those two classes of models is to compute the cosmic 
microwave background (CMB) anisotropies they predict and confront the results to 
observations. At present the anisotropies on large angular scales are known from 
the COBE measurements ||. On smaller scales, the present balloon and ground 
experiments [[| are still unprecise but the measurements should improve when the 
MAP @ and Planck [| satellite missions are launched. 

On the theoretical front, the CMB anisotropy predictions from the inflationary 
models (at least the simplest ones) are robust, and can now in many cases be 
obtained by plugging into the CMB FAST code || the initial conditions given by 
the particular model at hand. On the other hand the predictions from defect scenarios 
are still unclear, the reason being that the defects, which act as a continuous source 
of inhomogeneities, are difficult to model. Indeed the description of the formation, 
evolution and decay of defects requires heavy numerical simulations (see e.g. 
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and the CMB anisotropies their networks create seem to be very sensitive to their 
detailed structure. Recent results jl6]]|l7j] indicate that defect model predictions do 
not agree with the data. It is clear however that more work needs to be done before 
rejecting topological defects as a possible explanation for the CMB anisotropies as 
they may for example have a more complex internal structure (see e.g. JTsfl) than that 
assumed in those numerical models. 

Considering the complexity of modeling realistic topological defects, the semi- 
analytic approach initiated by Durrer and collaborators seems to be a fruitful 
compromise. The idea is to model the defect network by a stress-energy tensor which 
acts as a source to the linearized Einstein equations and induces inhomogeneities in 
the cosmic fluids. The ten components of this stress-energy tensor are then drastically 
constrained by a number of physical requirements. In this paper, we shall impose, as 
in (lijm, that the defect network 

(i) is statistically homogeneous and isotropic (hence obeys the copernican principle), 

(ii) is created at a phase transition in an up to then perfectly homogeneous 



and isotropic universe (hence obeys specific causality 21 and matching [|22f 
conditions), 

(iii) evolves deep in the radiation era and deep in the matter era in a way which is 
statistically independent of time (scaling requirement []l9| ) . 

We shall also make the further assumptions that the network 

(iv) satisfies the conservation equations, 

(v) is statistically coherent (see e.g. [p2| and below for a precise definition). 

Of all those requirements the latter two are stringent and probably not fulfilled by 
realistic defect networks (see e.g. @ f23| ) • 

Now, as for assumption (Q), it is a sensible one to make, as incoherent sources can in 
principle be described as a sum of coherent ones Q (see §|3.3|). As for assumption (Q, 
all authors who describe semi-analytically the defects as a fluid have up to now made 
it. All authors who consider global defects construct conserved stress energy tensors 
(see e.g. [pi] ). Those who study local strings may have to add an extra fluid so that 
the total stress energy tensor is conserved (see e.g. |l^][|4|). Only Contaldi et al. 
as far as we are aware, have questionned it. In order to take into account the network 
loss of energy through gravitational radiation or particle production, they have added 
an extra fluid which compensates for the non conservation of the defect stress-energy 
tensor (which comes from excising small loops from the numerical simulation in order 
for the network to scale). The problem with such a modelling (as Contaldi et al. 
themselves note) is, first, that the results are very sensitive to the equation of state of 
that extra fluid (which is left as a free parameter), and, second, that the energy lost 
by the network is transferred exclusively to that extra fluid and not to the background 
photons, baryons, Cold Dark Matter (CDM) and neutrinos. Therefore the modelling of 
the defect network energy loss is still in a primitive stage and needs to be improved by 
a careful analysis of the microphysics involved in the evolution of the network (which 
includes friction effects, particle and gravitational radiation production by decaying 
loops, as well as electromagnetic interactions with the plasma when the network carries 
currents). In this paper, in order to keep arbitrariness to a minimum, we choose to 
ignore the network energy losses. 
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Those five conditions being imposed, the stress-energy tensor of the defects will 
turn out to depend only on four free functions, two describing its scalar part and two 
describing respectively its vectorial and tensorial parts. Moreover, for a particular 
ansatz which garantees that the defects evolve from the two scaling regimes in such 
a way that the conservation equations are satisfied, those four functions will, as we 
shall see, depend only on one variable u and tend to constants when u — > 0, and to 
when u — > oo. 

We shall choose simple ansatze for those four functions (their precise shape can 
in principle be obtained by numerical simulations) and solve the well-known (see e.g. 
p5|) linearised Einstein equations which couple the defect network to the cosmic fluid 
inhomogeneities. Finally we shall compute numerically the CMB anisotropics those 
inhomogeneities produce. 

In agreement with [|l9| (as well as with |l6| [^0| ) we shall see that those anisotropies 
are very sensitive to the choices made for the free functions describing the defects. Our 
results therefore confirm previous work but the loss of scaling invariance during the 
transition is treated more carefully. They are also more complete as they include the 
vector and tensor contributions. 

The paper is organised as follows : in §|| we write down the linearised 
Einstein equations (in longitudinal gauge) describing the evolution of the cosmic 
fluid inhomogeneities driven by the defect stress-energy tensor (we suppose that the 
material content of the universe is a mixture of four perfect fluids : CDM, neutrinos, 
baryons and photons, the latter two fluids being tightly coupled before decoupling). 
In ^ we see in a detailed way how the five conditions given above constrain the defect 
stress-energy tensor. In we solve analytically the linearised Einstein equations 
deep in the radiation era, when all scales of interest today were larger than the 
Hubble radius, and we set the initial conditions on the cosmic fluid perturbations, 
following ||. In §| we give the correlation function of the CMB anisotropies in 
function of the cosmic fluid perturbations, assuming instantaneous decoupling (our 
computation is slightly different from the standard one — see e.g. p6fl). Finally 
in we compute numerically the multipole coefficients Cg of the CMB anisotropy 
correlation function using a perfect fluid code and compare our results to previous 
computations. We discuss the approximations made, to wit perfect fluid description, 
tight coupling and instantaneous decoupling hypotheses, by comparing our perfect 
fluid code to a Boltzmann one. We also show how to evaluate analytically the vector 
and tensor contributions to the Ci for small I. 



2. Einstein equations 

2.1. The background 

The universe at large appears to be remarkably homogeneous and isotropic and 
governed by the gravitational force created by its material content. We suppose that 
that was always so, and describe it in a first approximation by a Robertson- Walker 
geometry (that we shall for simplicity suppose spatially flat) whose time evolution 
satisfies Friedmann's equations. We suppose that its material content is a mixture 
of perfect fluids made of decoupled black body radiations of photons and massless 
neutrinos, of baryons with which the photons were coupled until their temperature 
dropped to w 3000 K, and of uncoupled "Cold Dark Matter" (CDM). We thus take 
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the line element to be : 

ds 2 = a 2 (r])(-dr) 2 + S lJ dx i dx j ), (1) 

where 77 is conformal time (we set the velocity of light equal to 1), where x 1 , i — 1,2,3 
are three cartesian coordinates and Sij is the Kronecker symbol, and where a(rj) is the 
scale factor. 

The conservation equations for the fluids read : 

p' n = -m( Pn + p n ), (2) 

where p n and P n are the energy density and pressure of the n-th fluid and where 
TL = a' I a is the comoving Hubble radius, with a prime denoting a derivative with 
respect to conformal time. The radiation fluids of photons and neutrinos both have a 
pressure equal to a third of their energy density, and the fluids of baryons and CDM 
have zero pressure (whence they have become non-relativistic). Hence the two former 
scale as a -4 and the two latter as a -3 . 
The Friedmann equation is : 

n 2 = Kjpt, (3) 

where k = 8ttQ is Einstein's constant and pt = Y^ Jn p n - 
Combining eqns. (||) and (||) yields : 

h 2 = ±(ti° r + xn° m ), (4) 

x z 

with 

n? = n° + 0°, (5) 
n° m = n° b + n° c , (6) 

x = a/ao (J) 
h = H/Ho, (8) 

where fl n = p n / Pt (with ^ n Q n = 1), where the index zero means the present time 
and where the indices 7, v, b, c refer respectively to the photons, neutrinos, baryons 
and CDM. 

The numerical values of are known from : 

(i) the value of the Hubble constant Hq = Tio/ao, which determines p® via eq. (||) 
(we shall take Hq — lOO/io km sec _1 Mpc~ 1 with ho = 0.5), 

(ii) the temperature of the microwave background (2.726 K), which determines p® 
via Stefan's law, 

(iii) electroweak theory, which determines the ratio p®/p® = 0.68f32 (assuming three 
families of massless, non-degenerate neutrinos), 

(iv) nucleosynthesis, which determines p° (we shall take f^/i 2 , = 0.0f25). 
All that yields : 

n° ~ 9.4 x fO^ 5 , (9) 
0° ~ 6.3 x f0~ 5 , (fO) 
n° b =5.0x IQ- 2 , (11) 

n° c ~ 0.95. (t 2) 
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The value x^ of x at decoupling is given by the ratio of the present temperature of the 
photon background and its temperature at decoupling : x^ « 0.9 x 10~ 3 . Equality of 
matter and radiation is defined by x cq = £!°/f^ and occurs before decoupling. 
Equation (El) can be easily solved to give x{rf) : 



Vo -V = 



1 - y/Sl° r +xW m . (13) 



m 



In the radiation era (x -C x e q), x oc rj, in the matter era [x S> x e q), x oc rj 2 and 
both eras H oc 1/r). 

Finally, when the scale factor behaves as a(rj) oc r\ q with q > when 77 — > (which 
is the case in the standard scenario we consider here) the particle horizon of the point 
x l is defined by r — rj (where r = \x l — x n \ is the comoving spatial distance from x l ) 
and is proportional to the comoving Hubble radius Hr 1 . 



2.2. Perturbation equations 

Following Bardeen |27J (see also p8[) we split the perturbations of the geometry and 
the matter variables into (spatial) scalar, vector and tensor components (for reviews of 
this formalism, see e.g. We shall work in the longitudinal (also called newtonian) 

gauge in which the line element of a perturbed Friedmann-Robertson- Walker space 
time reads : 

ds 2 = a 2 (r]){-(l + 2$)dr] 2 -2$ t dx\h]+[(l-2y)5 t3 + 2E lj ]dx l dx j }.(U) 

Here and in the following all barred spatial vectors (such as $i) are divergenceless 
(di& — 0), all barred spatial tensors (such as Eij) are divergenceless and traceless 
(diE 1 ^ = Q,El = 0) and all spatial indices are moved with the Kronecker symbol 
The six functions ($,\1/), $i and Eij are respectively the scalar, vector and tensor 
parts of the metric perturbations in longitudinal gauge. 

We shall treat the material content of the universe as a mixture of perfect fluids (see 
§|(] below for a discussion of this hypothesis) . The perturbations of the stress-energy 
tensor of the n-th perfect fluid in this perturbed universe are : 



<5T " = a 2 Pn (5 n + 2$) 

= a 2 [p n $ t - (v i>n + d l v n )( Pn + P„)] 
5T% = a 2 P n [2E l3 + 6ij(6 n - 2*)] 



(15) 
(16) 
(17) 



where 5 n = Sp n /p n is its density contrast and dtv n + Vi lTl its velocity perturbation. 

The stress-energy tensor of the topological defects is in itself a small perturbation 
(this is the so-called "stiff approximation", see e.g. 29|). We decompose its 
components as : 



K 



K 



e, 



p 



% + ;>m: + djUt 



ps 



-AU S 



d tJ n s 



(18) 
(19) 

(20) 



The ten source functions (p s , P s , v s , II s ), (w|,II|) and nf will be discussed later. 
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We shall write the evolution equations in Fourier space, the Fourier transform of 
any function f(x l ,rj) being defined (formally) as : 

f(k\v) = -?T¥ / d 3 * e - ife **" f{x\r,), (21) 



(27T) 

f(x\v) = / d 3 ke ikiX ' f{k\ri). (22) 

(The factor 7Yq insures that / and / have the same dimension.) 

As announced in the introduction, we assume that the defects interact only 
gravitationally with the rest of the matter. Therefore their stress-energy tensor is 
conserved, which yields for all (but the k % = 0) modes : 

= -Hp s -3HP S + kV s , (23) 

dr/ 

dV s ~ / - 2 \ 

— = - 2HV S + k l-P s + -TT s j , (24) 
= -2-Hvt + knl (25) 

where we have introduced k = \/kik l , where a barred vector is orthogonal to k l and 
where 

V s = kv s , (26) 
n s = k 2 fl s , (27) 

K = kill (28) 

After decoupling the fluids of neutrinos, photons, baryons and CDM are decoupled. 
Their stress-energy tensors are therefore separately conserved, which yields : 

(29) 
(30) 
(31) 
(32) 

U.7/ UV/ 

where we have introduced 

5l = L - 3 ( 1 + ^) 4>, (33) 
V Pn / 

V n = kv n . (34) 



tit 

d.77 


4 - 

= -w v 


dV v _ 
drj 


- fc 




dr] 


= 0, 


drj 


4 » 
= -kY, 


d% 
' dr] 


- fc 




dV^ 

dr] 


= 0, 


drj 


= kV b , 


dVb 
drj 


-nv b - 




dv b _ 
dx 


-Hv\, 


d^ 
drj 


= kV c , 


dVc 
drj 


-nv c - 


fed , 


d«l _ 
d?7 





Before decoupling the density contrasts and velocity perturbations of the neutrino and 
CDM energy fluids still obey eqns. ( p9| ) and (32) but the baryons and photons being 
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tightly coupled, eqns. (|3(]) and ([H]) must be replaced by (using also the fact that the 
number of baryons is conserved) 



di] 



Vh = Vy 



Wi,6 



l >7 



(I17 
df 7 

dr? 
d«j^ 

drj 



-HRVy 



(35) 
(36) 
(37) 



where we have introduced i? = (3ax)/(4 + 3aa;), a = (It is clear that imposing 

that the density contrasts and velocity perturbations of the baryon and photon fluids 
obey eqns. (|35-37) until decoupling and ( |30| , ^l| ) after decoupling, all variables being 
continuous at x = Xd, is a first approximation : see §[| below.) 

The Einstein equations give, after linearisation, the scalar and vector metric 
perturbations algebraically in function of the matter perturbations as 



* = 



-1 



dH/dq) 



$ = 



*i = 




1~L 

3(l + u n ) — V n 



k 



(38) 

(39) 
(40) 



(The other Einstein equations are redundant because of Bianchi's identities. The index 



n runs on the four background fluids, and we have set uib 



0J n 



0, LO v 



1 /3 and 



q — k/H.Q.) Finally the tensor metric perturbations are not determined algebraically : 
they solve the differential equation 



d 2 £, 
dry 2 



- 2H- 



dr) 



k 2 E i3 =nl%. 



(41) 



where a barred tensor is traceless and orthogonal to k % . 

The background quantity li.(r]) being known from eqns. ([|, |l^) and ( p9f{4l| ) give 
the evolution of the perturbations in function of n and the comoving wavenumber k 7 
once the source variables (constrained by eqns. (|23-25)), are known in function of r/ 
and k (Q and once the initial conditions are set (§ I[k 



3. The stress-energy tensor Mt/ of the sources 

3.1. Homogeneity and isotropy constraints 

The ten components of the stress-energy tensor Ml ,(?7, x l ) of the topological defects 
(/i, v = 0,1,2,3) are ten statistically spatially homogeneous and isotropic random 
fields. So are their ten Fourier transforms O^ij], k?) (which are complex but such 
that 8*^(77, fc 1 ) = Q^(r],—k 1 )). (We ignore the k % = mode which can be absorbed 
in the background.) 
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The statistical properties of those ten random fields will be described by their 
unequal time two-point correlators 



n. 



(42) 



where (...) means an ensemble average on a large number of realisations, and where 
the correlator C^po is a tensor which depends only on 77, 7/ and r % = x % — x n because 
of the spatial homogeneity of the distribution. The power spectra of the correlators 
Cfiupo- are defined as 



(43) 



where a hat denotes a Fourier transform. The power spectra are related to the 
correlators in Fourier space by 



(6* ( V , k^e pa (rj', k' 1 )) = 8(W ~ k H )^-P^ pa ( V , rf, k*). 



(44) 



The spatial isotropy of the distribution now forces the power spectra to be of the form 



-fbooo 


= A , 






(45) 


-Poooi 


= ihBi, 






(46) 


PoOij 


= CoSij + Cikikj, 






(47) 


PaiQj 


= D 5ij + D 2 k l k j , 






(48) 


PQijk 


= i [EikiSjt + Ei(kj5ik 


+ kkSij) 


+ E 3 kikjk k ] , 


(49) 


Pijhl 


= FoStjSki + F {5ik5ji + 


SuSjk) 4 


F 2 (kikj6ki + kkhSij) + 


(50) 




F 2 (k l k k S j i + k % ki8 jk 4 


- kjki8 ik 


f kjkkSu) + F^kikjkkh, 





where Aq, B\ etc are 14 real functions of r\, rj and the modulus k of the spatial vector 
k\ 

When Omj, is decomposed into its scalar, vector and tensor components according 
to eq. (|18||20|), then eqns. (g§) and @-|0|) yield 



(KM 

(pit) 

(P* S Ps) 

(v;v s ] 
(v;ps, 



(PIPs) 



(p;*s) 



A , 

-kB u 

Co + \k 2 C 2 , 

- k 2 C 2 , 
D + k 2 D 2 , 



= k 



E 1 + -E x 



— E% 

3 3 



= - k{2Ex + k 2 E 3 ), 



Fo + —Fq 



h 2 
3 



Fo H — Fo 

3 



Fo H — Fo 

3 



3F + Ak 2 F 2 + k 4 F A , 



^k 2 F 4 



(51) 
(52) 

(53) 

(54) 
(55) 

(56) 

(57) 

(58) 

(59) 
(60) 
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(CS->=A^i, (61) 

-kEkPa, (62) 

(%rr j ) = (F + k 2 F 2 )P ij , (63) 

(fty fti) = F (PikPji + PaPjk ~ Pijflw). (64) 



where (p* (77, k z )p s (r]', k n )) = S(k l - k n )(p* s p s ) etc, and where Py = - kikj/k 2 . All 
other correlators are zero so that the scalar, vector and tensor parts of fall into 
statistically independent sets. 

We also note that if (k 2 F 2 ) and (fc 4 Pi) are of higher order in k than F , then, for 
small k @|| 




3.12. Causality constraints 

Since the network of defects appeared at a definite time the distribution must be, for 
causality reasons, completely uncorrelated on scales larger than the particle horizon. 
Therefore, as stressed e.g. by Turok Q, the unequal time correlators are strictly zero 
outside the intersection of the past light-cones, that is : 

CW(?7,7/,r l ) =0 if r>n + rf. (66) 

Property ( |66| ) translates in Fourier space into the fact that the equal time power 
spectra are white noise on super horizon scales (that is for kr\ <C 1). Indeed, because 
the correlators (|6(]) have compact supports their Fourier transforms are C°° in k l . 
Therefore causality forces the fourteen functions B\ etc to be C°° in k 2 . 

Moreover, since within one horizon volume there are almost no defects, those 
fourteen functions must tend to zero on small scales, that is for krj 3> 1. 

3.3. Coherence hypothesis 

Any distribution of active sources such as topological defects must be homogeneous, 
isotropic and causal. The hypothesis of statistical coherence, that we shall now make, 
is stringent but motivated by the fact that any incoherent distribution can in principle 
be decomposed into a sum of coherent ones O] . 

By definition two statistically homogeneous random fields Si(r],k l ) and 5*2(77, A; 1 ) 
are statistically coherent if their correlators factorize, that is if 

(5:(r ? ,fc i )5 6 (r ? ' ) fc rt ))=5(fc < -fc' i K( J7 ,/ C >6(r ? ',A ; i ) (67) 

for a, b = 1, 2. For our purposes this is equivalent to saying that 

5 a (ry,F) = p a (ry, F)e(fc l ) (68) 

where p a (r], k l ) is a real function and where e{k l ) is a normalized complex random 
field such that (e*(fc J )e(fc' 1 )) = 8{k l - k' 1 ). 

Let us first combine coherence and isotropy, that is (]67|-p8|) and ([HJQ). The four 
scalar components p s , V s , P s , and n s of MI , form a statistically independent set. The 
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hypothesis ( pq ) of maximal coherence implies that 

p s ( V: k*) = p p { V ,k)e{k i ), (69) 

V s ( V ,k l ) = Pv (r ] ,k)e(k t ), (70) 

P s (n,k l ) =p P (r 1 ,k)e(k i ), (71) 

r(77,fc l ) = p 7r (r 1 ,k)e(k l ). (72) 

Hence, out of the fourteen functions A , B\ etc appearing in ( |5lH60| ) ten can be 
expressed in terms of the four real functions p a and the four functions (e.g. Do, E\, 
Fq and F2) that remain arbitrary. More precisely : 

Ao = P% (73) 

«--»(?), (»> 

Co =Pp\PP + ) , (75) 



(76) 

£1 = ^ (W + ^ , (77) 



/Py\ 2 _ Do 



<™> 

F, = -%(pr + ift,) + §, (80) 
Fo = (w + jP„) - Po, (81) 

The four vector components v i and tt^ of M „ form yet another statistically 
independent set. In order to obey ( |68| ) as well as (]6l|— 163|) they must fall into two 
independent subsets corresponding to two different polarisations. Hence we shall write 
them as 

%(T),k i )=Mv,k)e i {k i ), (83) 
r i (r,,k i )=p*(r 1 ,k)e i (k i ), (84) 
e^k 1 ) ^heiik^+fhiemik'), (85) 

where (e* A (k' t )eB{k t )) = 5AB$(k n — k l ) with A,B standing for / or m, and where l % 
and fh % form with k % jk an orthonormal basis (so that P^j = ZjZj + fhifhj). Hence two 
of the three previously free functions D , E\ and F can be expressed in terms of the 
two real functions p v and pjj as [see (61-53)] 

Do=P% (86) 

Ei = - ^p*, (87) 

F =pl-k 2 F 2 . (88) 
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Finally the two tensor components IL^ of M „ form a third statistically 
independent set. In order to obey (Q) they cannot be coherent but, on the contrary, 
independent of each other. Hence we shall write them as 

fh j (v,k i )=p T (v,k)e i j(k% (89) 
eijik 1 ) = %e+(fe i ) + 7%e x (fe i ), (90) 

where Uj - IJj — fhifhj and m,j = lifhj + Ijfhi are the two polarisation tensors, 
where (e^(fc' l )es(fc 1 )) = SABS(k n — k l ) with A,B standing for + or x, and where [see 
eqns. (p4) 1 

Fo=Pt- (91) 

In conclusion combining coherence and isotropy reduces the number of free 
functions from fourteen to seven : p p ,pv ,PPiP-k,Pv, Pit and pr- 

Let us add now the causality constraint which states that the fourteen functions 
Ao, B\ etc must be C°° in k 2 . It implies [see (|73"1)1 that p p , pv/k , pp and p^/k 2 are 
C°°. This in turn implies [see eqns. (f^fp)] that D /k 2 , E x jk 2 , F 2 /k 2 and F //c 4 
are analytic as well, and therefore, from (|86|), pc/fc and p^jk 2 are C°°. Finally, from 
©, p T is also C°°. 

An immediate consequence of this regularity requirement is that 

= V \ = 0{k% (92) 

(**%) -2p| = G(fc 4 ), (93) 

(ft" ff^> = 4p| = 0(fc 4 ), (94) 

so that, contrarily to the general case (^5|), the anisotropic stress correlators of coherent 
defects (when p n ^ p^ ^ p?) are not in a definite ratio for small k. 

A question that can be asked is how can a sum of coherent defects (with anisotropic 
stress correlators of order fc 4 ) leads to an incoherent defect (with anisotropic stress 
correlator of order k°). In a paper to be submitted Q, we show explicitely how the 
fc°-behaviour of incoherent defects can be obtained as a (infinite) sum of A: 4 -behaving 
coherent defects. We also show how the causality constraints fade away in the process 
and how the procedure introduces almost no extra arbitrariness. 

3.4- Scaling hypothesis 

We shall now make the other assumption that, as long as the evolution of the universe is 
described by a single scale (to wit the Hubble radius H^ 1 ), that is deep in the radiation 
era and deep in the matter era but not during the transition, the dimensionlcss 
quantities Q^/pt are scale invariant, that is such that 

C^u(v, V, r l ) = H 4 F^(r l H) (95) 

(using the fact that na 2 pt oc H 2 in both eras). This assumption is supported by 
a number of numerical simulations [0-[0(l£|[|3l| as well as qualitative arguments 
fl9f |32][ . In Fourier space this hypothesis translates as 

P^(V, V, k l ) = (2irfl 2 HF, iy {V/H). (96) 



This implies [see ( |45[]6l| )] that the functions A /Ti, D /Ti, F /Ti, F /Ti on one hand, 
D2H, F2H, F2H on another and finally are all functions of k/H only (when 
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rj = rj'). From eqns. (|5l[-p4[) one then obtains the behaviour of the equal time 

correlators (p*J s ), (V S *V S ), (P* P.), (P^s), (K* s ), (C^>> and (n^fty). 

Combining this scaling hypothesis with the one of coherence implies [see eqns. fl7q)- 
(p2|)] that all six p a /VH are functions of k/Ji only. 

Therefore, all in all, coherent defects, which scale during the radiation and matter 
dominated eras, which are causal and statistically homogeneous and isotropic, are 



described by the following random variables : 

P s =Vhf 1 (q/h,x)e(q t ), (97) 

P s =Vhf 2 (q/h,x)e(q i ), (98) 

V s = -S=f 3 ( q /h,x)e(q i ), (99) 

^ =^U(q/h,x)e(q l ), (100) 

=-^=f 5 (q/h,x)e i (q i ), (101) 

# « =^h(q/h,x)e l (q l ), (102) 

2 

% = ^Hi/ h ^)^M) ( 103 ) 

where we have reintroduced the dimensionless quantities x = a/ao, h = TL/TLq and 



q % = k l /Tto, and where the seven dimensionless functions f a (q/h,x) depend only on 
q/h deep in the radiation or matter eras, tend to constants on super-horizon scale 
(that is for kr\ oc q/h <C 1) and to zero on small scales. 



3. 5. Conservation laws and loss of scaling invariance during the radiation to matter 
dominated transition 

When the defects are described by eqns. (^7-103) the conservation equations (23-p5|) 
become 

1 flf„ r)f„ 3 9 

fa, (104) 
u 2 h, (105) 
(106) 

f a ux <± 

where we have introduced 

u = q/h, (107) 
1 0° 





o s df 3 

- 3uj)u— 

ou 


dx 


h§(34 


w)/ 3 


f 2 2 

= /a-3« 




- 6Ld)U— 

ou 


■ dh 
\- x— 

ox 


^~ 




= - 3/ 2 - 




- 6Ld)U— 

OU 


■ dh 
\- x— 

ox 


+ 


w)/ 5 


= w 2 /e, 



3 n° r + xty 



(108) 



Deep in the radiation or matter eras uj — > 1/3 or 0, and scaling solutions 
independent of x are readily found for ji , fa and fa in terms of the four free functions 
fx, fa and f$. During the transition from the radiation to matter dominated eras, 
scaling is lost. How can this loss of scaling be modelled, and how is the final result 
dependent on this modelisation ? 

Two different approaches have been followed : 
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(i) Durrer et al. pl| used a "sudden transition" approximation in which u> = 1/3 
until x — x cq and u> = afterwards. The problem here is that the functions f a 
cannot all be continuous at x = x cq and that energy momentum conservation is 
violated during the transitionf . 

(ii) Turok plL Hu et al. Cheung and Magueijo [[55) fix the time dependence of 
two source functions and determine the other two by integrating ( 104 - 106| ) . 

We follow a route similar to ([]) , which also amounts to making fairly artificial ansatze 
but which garantees that the conservation equations are exactly satisfied : we first 
impose 



fi = T{u), 
h = h(u)g{x), 
h = h{u)p{x), 
h = C(«). 



It then follows from (|0§^0§) that 



/s(«) 



3(3 + uj)g + 4xdg/dx 
u 2 T(u,x) 
3(3 + lu)p + Axdp/dx 



u Q(u,x) 



(109) 
(110) 
(111) 
(112) 



(113) 
(114) 



where Q(u,x) and V(u,x) are some functions, C°° in u 2 , that we impose to depend 
on u only 



g(u,x) = g(u) , v(u,x) = v( u ). 

This implies that g{x) and p(x) satisfy the following differential equations 

(3 + cj)g + \x^- = l-u 
6 dx 



(3 + uj)p- 



4 dp 
—x— 

3 da; 



(115) 



(116) 



(117) 



(up to irrelevant overall constants) with lo given by ( JlOq ) and g(x —* 0) = 1/5, 
p(x — > 0) = 3/10 (and g(x — * oo) = p(x — > oo) = 1/3). We hence arrive at the 
following expressions for the seven functions f a in terms of the four ar bitrary C °° 
functions ^(u), Q(u), V(u) and Q(u), with and p(ar) being given by ( 116 -117) : 



h = r, 



3lo)u 



dT 
du 



h 
h 



1 



u 2 g 



9, 



1 



(l + 3w) 



1 dJ 7 



2 4 
U U 



dQ 



(118) 
(119) 

(120) 
(121) 



(.9 ~ 1)--J~ ~ 6 ^ ~ : j 

i u era du 

f Note added in proof : In Durrer et al. |^], the loss of scaling is taken into account by interpolating 
the source functions computed respectively in the radiation and matter era. The fit is performed by 
using cither the function t e q/(t + ieq) or exp(—t/t eq ) and the results do not depend too much on this 
interpolation. However, energy conservation is still violated during the transition. 
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ft 



6 



1 



' 3 
h = 
fr = Q 



u 2 T>p, 

^ + (l + 3o;)| 



-Gg, 



1,, o \ ^ 



(122) 

(123) 
(124) 



In conclusion, we note that in our model both components of the defects stress- 
energy tensor are generically present, in contrary to e.g. the "pressure model" or the 
"anisotropic stress model" [b4J pa] . 



Initial conditions 



The set of equations (|29|-[l"l|) with the defects described by eqns. (97- 103| ) and 
124) will be solved numerically. However, deep in the radiation era when all scales of 
cosmological interest today were larger than the Hubble radius, an analytic solution 
exists, which is easily found when adding to the set (|2^-[ll]) the Einstein equations 
for scalar perturbations which we had not included because they were redundant (by 
virtue of the Bianchi identities). Those two extra equations can be written as : 



d* 

dx 



i 



V s 
2qh 



(125) 



d 2 ^ 

da; 2 



2qhx 2 
1 d* 

x dx 



[4(Q°K, + fi°t> 7 
9 



3x(n b v b + n° c v c )], 



4ft! 



x 2 I 2 m 



1 2 2 

3 qX 



(126) 



1 

q 2 x 
1 



rbN 




Deep in the radiation era, when xfl^ <C : hx ~ V^r ( an d x oc rj). The 
wavelengths we are interested in are then much bigger than the Hubble radius : 
q/h <C 1, that is q 2 x 2 <C O? (or u 2 <C 1). In that regime the conservation equations 
for the defects [eqns. (1 IS -124)] yield that the functions f a are constants such that 

(127) 
(128) 

(129) 



fi = -6f 2 +0(u 2 ), 
h = \h + 0(u 2 ), 



9 



u 2 h + 0(u 4 ). 



Let us first solve the equations for the scalar perturbations. The conservation 
equations for the fluids |32|, (35|p7|) impose that in the long wavelength limit all the 
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density perturbations o" n are constant, that is random variables independent of time. 

St « C\ (130) 

dl^C (131) 

^ « (132) 

^ « (133) 

where rs means deep in the radiation era and up to terms of order u 2 . Eq. ( |126| ) for 
W is therefore a closed equation in that regime. Its solution is 



*o + ^ + ^ 3/2 (/4 + / 2 )e(g I ) 



80° 



(134) 



(135) 



where "tp and are two integration constant random variables and where we have 
set f a = fi°. 3 ^ 4 lim u ^o fa- The perturbation $ then follows from eq. 



yo + ^ + ^m-ifMq 1 ) 



(136) 



The Euler equations Q29| , |32L |35| - |37| ) then give the velocity perturbations of the 
fluids 



'7 



<1 



'7 



1 



l 



2*o a: 



(137) 



80° 



^ c 



1 



(138) 



+ - a; 5 / 2 (2/ 2 -7/ 4 )e(^) 
63 



240°/ 



/o° 



4 7 



2*o he 



(139) 



+ ^ 5/2 (4/ 2 -5/ 4 )e(g l ) 



80" 



where VJ n , V^ 111 and V™ are constants of integration. The nine integration constants 

introduced (the four density perturbations 5^ , *o, *i and the thr ee initial velocity 
perturbations V™) are constrained by the remaining equations ( 125 ) and ( |38| ) 



*° = - ^ m +wn, 

v] n = - vl n , 



(140) 
(141) 



1G 









= v„ , 




- i in 




- V c 


v c 






X 


±i 


±i ~i in 


Vb = 





so that the system depends on seven integration constants as requested. 

Solving the equations for the vector perturbations, deep in the radiation era and 
for long wavelengths, is straightforward. The conservations equations ( p9| , |3^ , 35-|37|) 
yield 

(142) 

(143) 
(144) 



where v" n are three constant vectors orthogonal to k % , and eq. (|40|) gives 



x 



* 5/2 / 6 eV) 



(145) 



with $ 



= (2/ 9 2 )[3^c n + 4(^s; in 



I in 



Finally the tensor perturbations solve eq. 



)], and *! = 

0) [with n; 



, given by (pi] 



P'3 
-^0 



63 



(146) 



where Sg 3 and i?^ are two constant traceless tensors orthogonal to fc\ 

The solution we have obtained is valid deep in the radiation era, for wavelengths 
larger than the Hubble radius, but late enough so that the baryons and CDM are 
non-relativistic. At the epoch of the phase transition which gave rise to the defects 
all the matter was relativistic and the solution (given in fl22| ) is ( |1 30 )-|146[) where one 
sets fi° 



fl° = = and 0° = 0° = 1. As shown in |22j, the simple fact that the 



defects suddenly appear in an up to then perfectly homogeneous and isotropic universe 
in such a way that the standard general relativistic matching conditions are satisfied 
determines some of the constants of integration in terms of the f a . More precisely, 
two of the seven integration constants which appear in the scalar perturbation set, 
one of the three integration constant vectors for the vector set and the two integration 
constant tensors for the tensor set are determined by these matching conditions (see 
p2[ for details). The other constants must be fixed by the physics at the epochs when 
the baryons and CDM become non-relativistic. For example, if one assumes that the 
transition is sudden and occurs on a hypersurface of constant density then (see p2| ) 
VP must be continuous and therefore S b min = 0. 

It turns out however that almost all those initial conditions which fix the 
integration constants are "forgotten" at late times, that is long enough after the 
transitions. This is clearly seen for all the metric perturbations and the scalar velocity 



perturbations. Indeed, asymptotically [see eqns. (134-139 

2 



145 



146)] 



9^ /2 (/4 + / 2 )e(g l ), 
^ 3/2 (2/ 2 -7/ 4 ) e (<f), 



(147) 
(148) 
(149) 
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Eij 

Vc 



63 



2 1 „6/2 



45 

2 g 
63 



z°^(4/ 2 - 5/ 4 )e(<f ), 
x 5/2 (2/ 2 -7/ 4 )e(q 4 ). 



(150) 
(151) 
(152) 



In order now to obtain the asymptotic behaviours of the scalar density 
perturbations one must solve the equ ations at next order in q. This is easily done : 
injecting the asymptotic expression s (15l|T 52) for the scalar velocity perturbations 



into the conservation equations |32 , |35|-|37|) one gets 



K -> < 



;«5 



(153) 
(154) 



Finally the solutions (142-144) hold for all wavelengths so that the initial conditions 
for the vector velocity perturbations are not "forgotten" in the evolution (apart from 
v C i)- We shall take (rather arbitrarily) 



= 0. 



(155) 



In conclusion, we shall start the numerical integration of the set ( 29 - 4l[) with the 
initial conditions (147 T55| ), at x = x- m (in practice x- m = 10 -9 ), for various sets of the 
constants f%, fi, fe and (which will depend on our choices for the functions !F(u), 
G{u), V(u) and Q(u) which describe the defects [see eqns. flll8 - |l24 )]), and all relevant 
values of q (in practice 5 x 10~ 2 < q < 5 x 10 2 , which corresponds to wavelengths 
from 10 times to one thousandth of the horizon size today). 

An important point to note is that all scalar perturbations, collectively called 
S(r],k z ), are proportional to the same random variable e{k l ) such that (e* (k n )e(k 1 )) 
= 5{k n — k 1 ), and can be written as 



S(k l ) = S(fc)e(fc i 



(156) 



where S is a function of k only which is entirely determined once the two functions 
fi(u) and Q(u) describing the scalar components of the defects are known. As for the 
vector perturbations, collectively called Vi(k l ), they can be written as 



V. t (k l ) = V(k)e z (k l ) 



with [see Q] (e*(k H )e k (k 1 )) = 6(k H - k l )P lk . 
collectively called Tij(k l ), can be written as 

T ij (k i )=T(k)e ij (k i ) 

with [see ©] {e* 3 {k H )e kl {k 1 )) - 5(k H - k^P^ + P u P jk - P tJ P k i). The functions 

V(k) and T(k) are entirely determined once the functions 'P(u) and Q(u) describing 
the vector and tensor components of the defects is known. 



(157) 

Finally the tensor perturbations, 

(158) 
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5. The correlation function of the microwave background anisotropics 



The energy of a microwave background photon in the direction 7*, as measured by an 
observer O at Xq is : E~ = — {g^ lv K' l u v )\o where g Mt /|o is the metric ( pT| ) evaluated 
at Xq, K^\o (with g^K^K v = 0) is the wave vector of the photon at Xq and u^lo 
(with g^u^u" = —1) is the velocity vector of the observer O. E® depends on K°\q, 
on the direction of observation 7* = — (K' l /\K\)o and the position and velocity of 
the observer O. Assuming that the photon has followed a nul geodesic (no-collision 
hypothesis) since it has decoupled from the baryons on a surface of constant electronic 
density or, almost equivalently, constant photon density, it is a standard calculation 
(see e.g. [g6| for a detailed derivation) to relate E~ to E d = — {g^K^u^d, its energy 
at decoupling, which is independent of 7* when "measured" by an observer comoving 
with the baryons. The result is 



AT 



E° 



E° 



4> 



4> 



{div 



(159) 



4*> 



1 % {diV b - 



Vifi) 

. d 



'(,i 



where E° = E d Xd (xd = a(r]d)/a(riQ)), where the index means that the quantity is 
evaluated today at (770, Xq — 0), where the index d means that the quantity is evaluated 
at {j]d,x d = 7* (770 — rid)) and where the integrands in the integrals are evaluated on 
the photon trajectory that is at x^ — (77, 7* (770 — 77)). Assuming that the microwave 
background is and always was a perfect black body radiation, eq. (159) defines its 
temperature anisotropies. 



Eq. (159), where the reception event Xq is defined (rather arbitrarily) as belonging 
to a surface of constant photon density, defines AT/T\^ as a gauge invariant quantity. 
The first term is the sum of a monopole (—[<!> + -j<5 7 ]o) which does not depend on 7* 
and a dipole ([7* (diV + Vi)] ) which is a Doppler effect due to the peculiar velocity 
of the observer. In the following we shall ignore these contributions and define the 
temperature background anisotropies as 



62 

^7 



w 7 



(160) 



with 



w 7 



^7 



4> 



7 



Y v itb \d + / <J>-df7 



'(n 



($' + *') dry, 



'(,; 



>lfi 



% dTJ. 



(161) 

(162) 
(163) 



'(,/ 



The next step consists in expressing the perturbations in terms of their Fourier 
transforms [see eq. (|2l|)l and to expand e lkiX * as 



i^(2^ + l)jt{kx)P e (k.x/kx) 



(164) 
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where j$ is a Bessel function, Pg a Legendre polynomial, and where kx = \Jkik l \J XiX 1 
and k.x = fc^x*. Eqns. (161-163) then read (reintroducing the dimensionless wave 
vector q l = k 2 /H ) 



Q S ^ T = 



(27T) 



^ ^ i £ (2^ + 1) y ^dgdfy P*(rf,.^ Gf vr ' T (g < 



(165) 



where 5, V, T stand for either 5 1 , V or T, where il^- defines the direction of the vector 
q l in Fourier space, where n q = g/g and w here [integrating on x = a/ao rather than 
on 77 and recalling the notations ( |l 56^1581 ) ] 



Gf (<f ) = Gf{q)e{q% 
G^iq 1 ) = GY(qh%(q% 
Gj(<f ) ^ TO 7 V'e^), 



with 



TO =My d ) 



3e(y) 



v b \ d 



&xji{y)— (* + ¥) 
da; 



TO 
TO 



' i'<i 

■ ji(ya)vb\d 



r 1 d = 

- / dxji(y)—E, 
Jx d dx 



(166) 
(167) 
(168) 

(169) 



(170) 
(171) 



where y = kr = k(rjo — 77) = 2q(l — y/fl® + xO^J/Q^n [see eq. (ff3|)] and where the 
index "d" means that the quantity is evaluated at x = Xd- 

0§, 0Y, and 6l are random variables which fall into statistically independent 
sets so that the temperature anisotropy correlation function factorizes as 

= (e fl e ?2 ) - (e^efs) + (e£e£) + (e£e£> (172) 

where 71 .72 is the cosine of the angle between the two directions of observation 7J and 
1\- 



(173) 



Injecting eqns. (165 171) into eq. (|172j) one therefore gets 

(*i v > T ei v > T ) - 7^3 EH^ + + V 



q'dqdn ? Gf> V < T G S f V < T Pt(n q .7i)Pi(n q -j2)A 



S.V.T 



where 



.4 s =1, 



(174) 

-4 V - 727i'^j = 71-72 - K-7i)K-72), (175) 
A T = l l 2 lhh[{PikP 3 i + PaPjk - PijPki) (176) 
= 2^.72 - (7l g .7 1 )(r? g .7 2 )] 2 - [1 - (n 9 .7!) 2 ][l - (n q .j 2 ) 2 }. 

Using now the identity |37| 

47T 



/ 



dClq-P^fiq.yi) x Pt>(n q .^ 2 ) 



21 + 1 



Pl{ll-l2)5tt< 



(177) 
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one first obtains (see e.g. Q) 



c 



47T 

2 



q 2 dq|G 



S|2 



where Gf is the function of q given by eq. ( |169| ) . The other identity (j] 



(n q .71 ) P? (n g .71 ) (n q .72 ) P?' (n g .72 ) dSl^ 



(178) 
(179) 

(180) 



= 47T 
-47T 



+ 1 



(2£ + 3){2£+l) 

£ 



(2i+l)(2£- 1) 



1±1 

2<?+ l' 

2£+ 1 



^ + 2 

-Of'J-2 



21-3 



P?+i(7i-72) 
p?-l (71.72) 



together with standard recursion relations for the Legendre polynomials and Bessel 
functions |$7| then yield (see e.g. p6|) 



c 



'-_£(£ +1) / g 2 d<7 



Vd 



Vb\d + dx 



My) d 1 

7/ da; 



(181) 



(182) 



The computation of (©^©5,) can proceed along similar lines but is more involved. 



The result is (see e.g. pq]) 

(&l@l) = ^ £(2* + l)Cjp(7i.72) 



ith 



2)! 



7T (1-2)) 



The coefhcients C e ' ' , which will be computed numerically in the next section 
for various choices of the source functions describing the defects, are the quantities to 
be compared with observations, through the inverse formula 



q dq 



dx 



My) d g 

y 2 dx 



(183) 



(184) 



Cf + CY + Cj — 27r 
C(z) 



dzC{z)P e {z), 



7i-72, 

(%a&>. 



(185) 

(186) 
(187) 



(Observationally, the ensemble average (...) is replaced by an average on different 
regions of the celestial sphere. For large angles there are few such independent regions 
so that the "cosmic variance" problem arises ; see e.g. pp| .) 

f We thank Luc Blanchet for showing us how to demonstrate it by means of the STF formalism (see 
e.g. M) 
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6. Numerical results 



Two of us (JPU and AR) have developped a numerical code to solve the equations for 
the perturbation s (p9[j4l] ) and compute the coefficients Cg in function of the multipoles 
i [eqns. ( |I79|J182|JL8^ ] when the defects are described by eqns. ^-|0§ and (|118Hl24| ), 



and the initial conditions are given by eqns. (147-155) (the computations start at 
x hl = 10~ 9 , q e [2 x 1(T 2 - 2 x 10+ 3 ] and I G [2 - 1500]). This code incorporates Silk 
damping (that is the decay of the Cg at large £ due to the fact that a fluid description 
no longer holds on small scales) by a simple multiplication of G s in Q169| ) by the 
decaying exponential (see 0|-E3]) 



-2fc ni, 



(188) 



where x s = OmU^ttbho)- 1 / 2 ^^ @|. 

As for the fact that decoupling is not instantaneous, we take it into account 
by convolving eqns. (169-171) with the visibility function (see Eg] for a detailed 

i 2 2 2 

derivation; a simple multiplication by e ' ) « ! <= cr with a — 0.03 for standard 
recombination p2|, is too rough an approximation) . 

Instead of eq. ( |188| ) we could have used the better fit proposed by Hu and Su giya ma 
|E3 but it is an analytical fit deduced from a Boltzmann code whereas eq. ( |188| ) is 
obtained from physical arguments flij]| . 

In this section we first test our code and discuss the approximations made. We then 
give a few examples of solutions and compare the results with previous calculations. 
(Note that for a better visibility, the results are not normalized to the COBE data at 
I = 10 as is usual 

6.1. No sources and adiabatic conditions : comparison with a Boltzmann code 

The code which we have implemented can of course also be used to study scenarios 
where no defects are present and where adiabatic initial conditions (as imposed 
by e.g. the standard inflationary scenario) are imposed. One has just to set 
J-=G = V= Q = and replace the initial conditions ( [147 -155) by 



^ = S\ = ~$ = = - 6*0 (189) 



3 3 

1 qx 



V v = Vy = V b = V c = - 77^=*o (190) 



# = $ =* (191) 

where ^o, in the standard chaotic inflationary model, is fixed by the quantization of 
the inflaton so that (see e.g. p6[ ) ^ cx q~ 3 ^ 2 . (We restrict ourselves to the scalar 
perturbations as the vector and tensor contributions are negligible in the standard 
inflationary scenario — see e.g. |47j].) 

The Cf coefficients as computed with our fluid code, with the Boltzmann code 
we have developped |50|, and with CMBFAST jly] are compared on fig. [j] (top). One 
first notes the excellent agreement between our Boltzmann code with CMBFAST. One 
also notes that if the position of the Doppler peaks, as given by our fluid code and 
the Boltzmann ones, coincide, their heights differ by about 15 %, which we consider 
enough for our present purposes considering that the physics in our fluid code (as well 
as in e.g. Seljak's in p2|) and a Boltzmann code differs on four points : 
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We use a perfect fluid approximation to describe the various matter and radiation 
components of the universe, 

We make a tight-coupling approximation, that is we assume that the baryons and 
photons have exactly the same velocities until decoupling, 
We suppose instantaneous decoupling, 

We correct for those assumptions by an adhoc damping of the scalar perturbations 



on small scales [see eq. (188)] and we take into account the finite width of the 
last scattering surface. 

whereas a description based on Boltzmann's equation is used in our Boltzmann code 
and in CMBFAST. We find that the difference in the height of the peaks is due to 
these approximations (see [|oJ for a detailed discussion). 

In conclusion a perfect fluid code is useful to find the positions of the peaks 
predicted by various types of topological defects but can be inaccurate when it comes 
to their height. In the following we shall use a perfect fluid code and therefore comment 
only on the position of the peaks. In |50) we shall use our Boltzmann code and shall 
therefore be able to compare not only the positions but also the heights of the peaks 
predicted by various defects as well as with those predicted by inflationary scenarios, 
the latter being also computed using a Boltzmann code. 

6.2. The scalar contribution of coherent defects to the CMB anisotropics 



The scalar component Cf of the Ce coefficients is given by ( |179[ ) and is entirely 
determined once the two functions T(u) and Q(u) (with u = q/h = k/Ti) describing 
the scalar component of the stress-energy tensor of the defects are known. 

First fig. U shows, on the particular case Q{u) = and T(u) cx exp(— u 2 ) (fig. ||, 
bottom), that, contrarily to what happens in inflationary scenarios (fig. ^, top), 
the integrated Sachs- Wolfe (ISW) contribution to Cf (that is the integral term in 
eq. ( |169| )) dominates at low I. This is in agreement with intuition : in inflationary 
scenarios the large scale anisotropies are built in the initial conditions, whereas in 
defect scenarios they are due to line-of-sight gravitational effects. We are aware that 
the inverse Fourier transform of a gaussian does not meet condition ( p6| ) but we note 
that the Fourier transform of any compact supported C°° function decreases faster 
than any power law. 

Fig. |^ illustrates the case when = and T = cxp(— u 2 /L 2 ) for various values of 
L. The most striking feature of these curves is that they do not exhibit a plateau at 
low t when 1^1. As I decreases the first "peak" is shifted to lower multipoles and 
tends to the position predicted by the standard inflationary scenario (see fig. 0, top). 
This can be explained by the fact that the defects decay on larger scales. 

Fig. |i| illustrates the case when both T and Q are proportional to exp(— u 2 ), with 
Q/J- = r. First, as expected, when r is small one recovers fig. ||, bottom or fig. || with 
L = 1. Second, when Q ^ 0, the dominant component to the stress-energy tensor of 
the defects is their anisotropic stress II s [see eqns. ( 118421 and p7|)]. The first Doppler 
peak is then washed out. This result compares to [|17| where incoherent sources are 
studied, for which the anisotropic stress dominates too. An interesting particular case 



is when a ~ 0.3 which minimizes the anisotropic stress contribution [sec eq. (121)]. 
The first Doppler peak is then higher than the plateau. 

Finally, as a caveat not to comment too seriously on the height of the Doppler 
peaks when a perfect fluid code is used, fig. [I], bottom compares Cf as given by the 
perfect fluid code used here and the Boltzmann code used in Q for Q{u) = and 



23 



5 - 



4 - 



3 - 



2 - 




1 ' ' — ' ' ' — ' ' ' — 

10 100 1000 

I 



1.4 - 
1.2 - 




10 100 1000 

I 



Figure 1. Comparison between our fluid code (solid line), our Boltzmann code 
(dashed line) and CMBFAST (short-dashed line) in the case of inflation (top) ; 
Comparison between our fluid code (solid line) and our Boltzmann code (dashed 
line) in the case of topological defects (bottom) with T = exp(— u 2 ) and Q = 0. 
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Figure 2. Decomposition of the different contributions of the Cf in the case of 
inflation (top) and topological defects (bottom) with T = exp(— u 2 ) and Q = 0. 
In each case, we have represented the Sachs- Wolfe, Doppler, Integrated Sashs- Wolfe 
contribution, as well as the total contribution (dashed, short-dashed, dotted and solid 
lines respectively). 
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Figure 4. The effect of the function Q = exp(— u 2 ). For that purpose we have varied 
the ratio r = Q jT with r = 0, 0.1, 0.3 (dotted, dashed and solid line respectively). 
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J~{u) cx exp(— u 2 ). The much bigger discrepancy in the height of the peaks in that 
case as compared to the case of inflation is due to the fact that eq. ( [L88D for Silk 
damping is valid only if the metric perturbations in the evolution equations for (5 7 
can be aproximated by constant (see e.g. eq. (98) of Hu and White [|6|). It can be 
shown that this is approximatively true in the case of inflationnary scenarios in a flat 
universe without cosmological constant, but here it is not true, as it can be checked by 
looking at the ISW contribution. Of course, it is also possible to do some a posteriori 
corrections by comparing fluid and Boltzman codes and by fitting ad hoc damping 
terms in order to find a good agreement between the two, see Hu and Sugiyama p6fl , 
but we did not find them necessary since we simply wanted to present a semi-analytic 
formalism. Let us note finally that in the case of the pressure model p^||^5[ , the 
discrepancy between this (simple) fluid code and an accurate Boltzmann code is less 
serious, probably because the ISW contribution is smaller. 



6.3. The vector and tensor contributions 

The vector and tensor components of the CV coefficients [given by eqns. 
are entirely determined once the functions V(u) and Q{u) describing the vector and 
tensor components of the stress-energy momentum of the defects are known. Since we 
supposed v\ = they both reduce to their integrated Sachs- Wolfe contributions. 



0.001 - 



0.0001 



Figure 5. The vector (dashed line) and the tensor (solid line) contributions to the 
CMB temperature anisotropics in a model where V = Q = exp(— u 2 ) 



6.4- Analytic estimate of the C^ ,T at small £ 

As we have seen, only the integrated Sachs- Wolfe effect contributes to the vector 
and tensor parts of the Cg. In this section, we estimate them at small I. We first 
assume that the last scattering surface is deep enough in the matter dominated era to 



approximate Tt and u by : TL = 2/ij and u — kr//2. The equations (|40|) and 
$i and E'fci can then be rewritten [using eqns (|25|,101-103) and (157-153)] as 



= / 4 = 
$ + -$ 
V 

~/i 4 -' 
E + -E 
V 



= « —J6{U) 



k l E = K 



V2 
krj 3 / 2 
2V2 



Mu) 



the solutions of which are given by 



5. a 
$ = — r 

u 4 



jU-3/2 /•« 

' « U / 2 /6(«)^, 



A = a h p x 

u u 2u 
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(192) 
(193) 

(194) 
(195) 



- ,3 /2 jiMmW-jiWmM /7(u)du _ 



We now make the ansatz that the defects disappear as soon as they enter the 
horizon i.e. that fe(v) oc /7(f) oc Y(v — 1) where Y is the Heavyside function. 

Computing (194-195), their time derivative and inserting the result in 
eqns. (182,183), we obtain that 



c[ oc e(t + 1) 



y 2 + 1 



(1+2)! ( A fM 1 
(*-2)l./ y 4 y + 1 



dy, 



(196) 
(197) 



with A = r/o/rid — 1. The main contribution to the two integrals comes from the 
points y ~ I (since the Bessel functions are peaked around that point). In a matter 
dominated universe, we have that A = \J clq / ad — 1 — 32 and for I < 30 we will assume 
that A ~ 00. With that approximation, we are led to compute the integrals 



B = 



3_M 1 
yP" 1 y + 1 



dy 



-Kp<2£ + 2, 



(198) 



with p — 5 and = 3 for the tensor and vector modes respectively. Using the relation 



jKv) 



dy 



o 

tt r(i+ P ) r(£ + i/2- P /2) 

2 2 +p r 2 (l+p/2)r(f + 3/2 + p/2) ' 



(199) 



-Kp<2£+1, 



it can be shown that 

B e p = 4(1 - 7$ 



with 



1 



3.2 11 (I - 5/2) : 
9?r 1 



°^3< 25 (£ _ 3/2) 



(200) 

(201) 
(202) 
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Thus, the i dependance of C e ' at small £ is given by 



£(£ + l)Cj oc 



£(£ + l)Cj cx 



1) 



[i+2){i-iy 
i(t + 1) 

(£ + ?>)(£ -2) ' 



3 < £ < 20. 



(203) 
(204) 



This solution is compared to our numerical results on fig. || where Cj /Cj and its 
analytic estimate are plotted (when f$ — f-j). The agreement is excellent for 3 < I < 20 
as expected. 

This estimation only assumes that we are in a matter dominated universe, that 
Vo/Vd — 1 < 20 and that f$ oc fj oc Y. The two first approximations are well satisfied 
for small multipoles. The ansatz concerning the mod elisation of the function s /e and 
fi will affect only the numerical pre- factors in eqns. (196-197) and (|203 -204), as long 
as the scale of decay of the defects is comparable with the horizon. 



o' o 




Figure 6. Comparison between the vector and tensor part to the CMB anisotropics 
(solid and dotted lines respectively) wit h their analytical analytic estimates (+-es 

4. It can be seen that the agreement is good 
> 30 for the tensor part. 



and x-es respectively) computed in 
up to a multipolc of I ~ 20 for the vector part and 



6.5. Conclusions 

In this article, we have studied the signature of coherent scaling defects on the 
cosmic microwave background. For that purpose, we have developped a semi-analytic 
formalism to describe the topological defect network as well as two numerical codes 
(based respectively on a fluid and a Boltzmann approach) to compute the scalar, 
vector and tensor components of the coefficients Ci of the CMB anisotropy correlation 
function. 
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Concerning the scalar component, and contrarily to what happens in standard 
inflationary scenarii, there is an important contribution from the integrated Sachs- 
Wolfe term, which does not in general build up as a plateau. In the generic case, this 
term dominates at low I so that the power on large scales is of the same order or even 
higher than on small scales. Such characteristics do not fit the present observations 

However, one can improve the model in various ways. First one can find special 
combinations of the arbitrary functions which yield results that are more compatible 
with the data (see fig. [| and p0[). But of course such a fine tuning needs to be 
justified. Second one can play with the cosmological parameters. Third one can 
extend our approach to more realistic defects (incoherent defects, loss of scaling ...). 
Finally one can consider defects that are produced at the end of an inflationary period 
p3| (and not, as in this paper, in an up to then perfectly homogeneous and isotropic 
universe) . Various extensions of the standard model of particle physics predict such 
a possibility (see e.g. f48|| ). We plan to explore those various improvements in future 
work HJ. 
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